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Abstract 



^\i . As the second stage of the project multi-indexed orthogonal polynomials, we present, 

^ I in the framework of 'discrete quantum mechanics' with real shifts in one dimension, the 

muhi-indexed (g-)Racah polynomials. They are obtained from the (g-)Racah polyno- 
00 ■ mials by multiple application of the discrete analogue of the Darboux transformations 

•^ i or the Crum-Krein-Adler deletion of 'virtual state' vectors, in a similar way to the 

cn \ multi-indexed Laguerre and Jacobi polynomials reported earlier. The virtual state 

vectors are the 'solutions' of the matrix Schrodinger equation with negative 'eigenval- 
ues', except for one of the two boundary points. 



X : 1 Introduction 

This is a second report of the project multi-indexed orthogonal polynomials . Following the 
examples of multi-indexed Laguerre and Jacobi polynomials [1], multi-indexed (g-)Racah 
polynomials are constructed in the framework of discrete quantum mechanics with real shifts 
[2]. It should be emphasised that the original (q'-)Racah polynomials are the most generic 
members of the Askey scheme of hypergeometric orthogonal polynomials with purely discrete 
orthogonality measures [21 HI O [6] . They are also called orthogonal polynomials of a discrete 
variable [7]. These new multi- indexed orthogonal polynomials are specified by a set of indices 
V = {di, . . . , c^m} consisting of distinct natural numbers dj G N, on top of n, which counts 
the nodes as in the ordinary orthogonal polynomials. The simplest examples, V = {£}, i > 1, 
{Pi^n{x)} are also called exceptional orthogonal polynomials [8]-[27]. They are obtained as the 



main part of the eigenfunctions (vectors) of various exactly solvable Schrodinger equations 
in one dimensional quantum mechanics and their 'discrete' generahsations, in which the 
corresponding Schrodinger equations are second order difference equations [21 [281 129] • They 
form a complete set of orthogonal polynomials, although they start at a certain positive 
degree (£ > 1) rather than a degree zero constant term. The latter situation is essential for 
avoiding the constraints of Bochner's theorem p^. The exceptional Laguerre polynomials 
with two extra indices T> = {di, ^2} were introduced in [31]. We are quite sure that these new 
orthogonal polynomials will find plenty of novel applications in various branches of science 
and technology as other orthogonal polynomials. One obvious application is the birth and 
death processes [32]. These new orthogonal polynomials provide huge stocks of exactly 
solvable birth and death processes [33]. The transition probabilities are given explicitly, not 
in a general spectral representation form of Karlin-McGregor [34J. An interesting possible 
application is to one-dimensional spin systems and quantum information theory [35j. 

The basic logic for constructing multi-indexed orthogonal polynomials is essentially the 
same for the ordinary Schrodinger equations, i.e. those for the Laguerre and Jacobi polyno- 
mials and for the difference Schrodinger equations with real shifts, i.e. the (g-)Racah poly- 
nomials, etc. The main ingredients are the factorised Hamiltonians, the Crum-Krein-Adler 
formulas [361 1371 138] for deletion of eigenstates, that is the multiple Darboux transformations 
[39] and the virtual states solutions [1] which are generated by twisting the discrete symme- 
tries of the original Hamiltonians. Most of these methods for discrete Schrodinger equations 
had been developed [21 [261 [23 [23 [iQl [HI [12] and they were used for the exceptional (g-)Racah 
polynomials [23] • The concept of virtual state 'solutions' requires special explanation in the 
present case. In the ordinary quantum mechanics cases, the virtual state solutions are the 
solutions of the Schrodinger equation but they do not belong to the Hilbert space of square 
integrable functions due to the twisted boundary conditions. In the present case, the Hamil- 
tonians are finite-dimensional real symmetric tri-diagonal matrices. Therefore the eigenvalue 
equations for a given Hamiltonian matrix cannot have any extra solutions other than the 
genuine eigenvectors. Thus we will use the term virtual state vectors. As will be shown 
in the text, virtual state vectors are the 'solutions' of the eigenvalue problem for a virtual 
Hamiltonian Ti', except for one of the boundaries, x = Xmax fl2.50l) . The virtual Hamiltonians 
are obtained from the original Hamiltonian by twisting the discrete symmetry and they are 
linearly related to the original Hamiltonian f l2.46p . Thus the virtual state vectors 'satisfy' 



the eigenvalue equation for the original Hamiltonian, except for one of the two boundaries. 
The polynomial part of the virtual state vectors had been used for the exceptional (g-)Racah 
polynomials. One distinctive feature of virtual states deletion in discrete quantum mechanics 
with real shifts is that the size of the Hamiltonian matrix (xmax) remains the same. This 
is in marked contrast with the eigenstates deletion (Christoffel transformations [3lll2]), in 
which case the size decreases by the number of deleted eigenstates. 

This paper is organised as follows. In section two, the main ingredients of the theory, the 
difference Schrodinger equation for the (g-)Racah system, the polynomial eigenvectors and 
virtual state vectors, are introduced. Starting from the general setting of discrete quantum 



mechanics with real shifts in § 12.11 the basic properties of the (g-)Racah systems are recapit- 
ulated in § 12. 2[ Based on the twisting (symmetry), the virtual state vectors are introduced 
in § 12.31 Section three is the main part of the paper. The basic logic of virtual states deletion 
in discrete quantum mechanics with real shifts in general is outlined in § 13.11 The explicit 
forms of multi-indexed (g-)Racah polynomials are provided in § 13.21 The final section is for 
a summary and comments. 

2 Original System 
2.1 General formulation 

Let us recapitulate the discrete quantum mechanics with real shifts developed in [2]. We 
restrict ourselves to the finite dimensional matrix case, Xmax = ^• 

The Hamiltonian "H = {'Hx,y) is an irreducible (that is, not the direct sum of two or 
more such matrices) tri-diagonal real symmetric (Jacobi) matrix and its rows and columns 
are indexed by non-negative integers x and y, x,y = 0, 1, . . . , Xmax- By adding a scalar 
matrix to the Hamiltonian, the lowest eigenvalue is assumed to be zero. This makes the 
Hamiltonian positive semi-definite. By a similarity transformation in terms of a diagonal 
matrix of ±1 entries only, the eigenvector corresponding to the zero eigenvalue can be made 
to have definite sign, i.e. all the components are positive or negative. Then the Hamiltonian 
Ti has the following form 



n.,y = -^/B{x)D{x+l)6,+i,y - ^ B{x - l)D{x) d^.i^y + {B{x) + D{x))5x^y, (2.1) 
in which the potential functions B{x) and D{x) are real and positive but vanish at the 



boundary: 

B{x) > (X = 0,l,...,Xmax- 1), -B(Xmax) = 0, 

D{x)>0 (x = l,2,...,x^ax), D{0) = 0. (2.2) 

The Schrodinger equation is the eigenvalue problem for the hermitian matrix Ti, 

'H(t)n{x) = £n(t>n{x) (n = 0, 1, . . . , n^ax) , = ^q < ^1 < " " " < ^n^ax) (2-3) 

where the eigenvector is (f)n = (0n(a^))x=o,i,...,xmax ^^^ ^max = N. Reflecting the positive 
semi-definiteness and based on the boundary conditions (12 .2^ . the Hamiltonian (12.11) can be 
expressed in a factorised form: 

Ax,y = \/B{x) 6x,y - ^/D{X+ l)6x+l,y, {A^)x,y = \/B{x) 5x^y - a/D(x) 6x-l,y (2.5) 

Here A {A^) is an upper (lower) triangular matrix with the diagonal and the super(sub)- 
diagonal entries only. The zero mode equation, A4>q = 0, is 



^B{x) Mx) - Vd{x + 1) M^ + 1) = 


= {x-- 


= 0,1,. 


■ ) 2^max 


-1), 


(2.6) 


\/5(a:i„ax) 0o(a;max) = 0, 


(2.7) 



and the second equation is trivially satisfled by the boundary condition -B(xmax) = 0. The 
groundstate eigenvector is easily obtained: 



(f)o{x) 



X~l J~y/ \ 

\ W tk In (^ = 0,l,...,x^ax), (2.8) 



with the normalisation 0o(O) = 1 (convention: nfc=n * ~ -'-)■ Needless to say it is positive 
for X = 0, 1, . . . , Xmax- For the explicit examples treated in [2], <Pq{x) can be analytically 
continued to the entire complex x-plane as a meromorphic function and it vanishes on the 
integer points outside the boundary; (pl{x) = (x G Z\{0, 1, . . . ,Xmax})- The eigenvectors 
are mutually orthogonal: 






= ^ <Pn{x)(j)m{x) = -^Snm (ri, m = 0, 1, . . . , n^ax) • (2.9) 

1=0 " 



For simplicity in notation, we write Ti, A and A^ as follows: 

e^^ = ((e±^),J (x,?/ = 0,l,...,x^ax), {e^%,y = 5,±i,y, {ey = e'^, (2.10) 



U = -y/B{x) e'^VD(x) - y/D{x) e-'^^B{x) + B{x) + D{x) 



= -^B{x)D{x + 1) e^ - ^B{x - l)D{x) e~^ + B{x) + D{x), (2.11) 

A = ^B{x) - e^y^D{x), A^ = ^B{x) - ^ D{x) e"^. (2.12) 

For the Schrodinger equation (I2.22p . it is sufficient that the functions B{x), D{x) and 
0„(x) are defined only for the integer grid, x = 0,1, . . . , Xmax- In this paper we consider the 
case that the potential functions B{x) and D{x) are rational functions of x or (/^ (0 < g < 1). 
So they are defined for any x G C (except for the zeros of their denominators) , see the explicit 
forms (I2.20l) - (l2.2ip . Also we consider the eigenvectors in a factorised form: 

M^) = M^)Pnix), Pnix) =Pn{vix)). (2.13) 

Here Pn{v) is a polynomial of degree n in 77 and the sinusoidal coordinate ri{x) is one of the 
following [2J; //(x) = x, e'x{x + d),l — q^, q~^ — 1, e'{q^^ — 1)(1 — dq^), (e' = ±1). Since Pn 
is a polynomial, Pn{x) is defined for any x G C The Schrodinger equation (I2.22p gives a 
square root free difference equation for the polynomial eigenvector P„(x), 

S(x)(P„(x)-P„(x + l))+D(x)(/'„(x)-P„(x-l)) =^,P„(x) (VxgC). (2.14) 

2.2 Original (g'-)Racah system 

Let us consider the Racah (R) and the g-Racah (gR) cases. We follow the notation of [2]. 
Although there are four possible parameter choices indexed by (e, e') = (±1, ±1) in general, 
as explained in detail in § V.A.I and § V.A.5 of [2j, we restrict ourselves to the (e, e') = (1, 1) 
case for simplicity of presentation. The set of parameters A, which is different from the 
standard one (a, (5, 7, 5) 0, its shift 5 and k are 

R: X ={a,h,c,d), 5 = {1,1,1,1), «; = 1, (2.15) 

qR: q^ = {a,b,c,d), 5 = (1,1, 1,1), n = q-\ < q < 1, (2.16) 

where q^ stands for g'-'^i'^^S'---) = (g'^\ q^^, ■ ■ •)• We introduce a new parameter d defined by 

TdcfJa + 6 + c-(i-l:R (0 ^7\ 

'^~ \ abcd-\-' :gR ' ^^'^'^ 

We adopt the following choice of the parameter ranges: 

R: a = -N, < d < a + b, < c < 1 + d, (2.18) 



qR : a = q , < ab < d < 1, qd < c < 1, 



(2.19) 



and Xmax = "^max = N. They are sufficient for tlie positivity of B{x; A) and D{x; A) below. 
Here are tlie fundamental data [2]: 

{x + a){x + b){x + c)(x + d) 

B{x;X) 



:R 



■.qR 



D{x;X) 



{2x + d){2x + 1 + d) 
(1 — aq^){l — bq'-^){l — cg^')(l — dq^) 

(l-rfQ2x)(l-dg2x+l) 

{x + d — a){x + d — b){x + d — c)x 

{2x - I + d){2x + d) 
~ (1 - a-Mg^)(l - rMg^)(l - c-irfg^)(l - g^') 



:R 



-rf 



'H{X)(j)n{x] A) = £:n(A)0„(a;; A) (x = 0, 1, . . . , Xmax; n = 0, 1, . . . , ri; 
0„(x; A) = 0o(x; A)P„(x; A), 



:gR 



Sn{X) = 

V?(a;; A) 



n(n + d) 



:R 



(g""- 1)(1 -rfg") : gR 
2x + d+l 



, ri{x] A) 



x(x + d) : R 

(g-^-l)(l-dg^) : gR 



rf+1 



1-rfg 
P„(x;A) = P„(77(x;A);A) 



R 



: gR 



4-^3 



— n, 77, + rf, —X, X + d 



4</'3 



Ij :R 
g ; g ) : gR 



a, b, c 
'g~", dg", g~^', dg"^ 

a, 6, c 

J i?„(?7(a:; A); a — 1, d — a, c — 1, 0? — c) :R 

\ Rn[l + d + 7]{x; X); aq"^ , da^^ , cq^^ , dc^^\q) : gR 
[a, 6, c, d)x 2x + d 



f>o{x-xy 



dn{Xf 



{1 + d — a,l + d — b,l + d — c, l)x d 
{a,b,c,d;q)x 1 - dq"^"" 

{a~^dq,b^^dq,c~^dq,q;q)xd'^ 1 — (i 

{a,b,c,d)n 2n + d 

[1 + d — a,l + d — b,l + d — c, 1)„ d 
(-1)^(1 + d-a,l + d-b,l + d-c)N 

V" 

{d+l)N{d+l)2N 

{a,b,c,d;q)n 1 - dg^" 



R 
gR 



(2.20) 

(2.21) 

(2.22) 
(2.23) 

(2.24) 
(2.25) 

(2.26) 

, (2.27) 
(2.28) 



R 



{a~^dq,b~^dq,c~^dq,q]q)nd'^ 1 — d 

■l)^(a-idg, b-^dq, c-^dq ; g)iv (i^g^^(^+^) 



X- 



gR 



(2.29) 



((ig;g)Ar((ig;g)2iv 
Here /?„(• ■ ■ ) in (12.271) are the standard notation of the (g-)Racah polynomial in [5]. It should 
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be emphasised that the quantities B{x;X), D{x;X), £^„(A), -P„(x; A), (f)o{x; X)'^, c^n(A)^ are 
formally symmetric under the permutation of {a,b,c), although their ranges are restricted 
as above by (l2J8D - (l219l) . 

Here is a remark on the polynomial -P„(x; A), which is in fact a polynomial in the sinu- 
soidal coordinate 77(2;; A) fl2.24p . The sinusoidal coordinate has a special dynamical meaning 
[21, [28], Us]. The Heisenberg operator solution for //(x; A) can be expressed in a closed form. 
This means that its time evolution is a sinusoidal motion. Let R be the ring of polynomials 
in X (the Racah case) or the ring of Laurent polynomials in q^ (the g-Racah case). Let us 
introduce an automorphism X in i? by 

I(x) = -x-d : R, X(g^) = q-'^d-^ : qR. (2.30) 

Obviously it is an involution X^ = id. The following remark is important. 

Remark: If a (Laurent) polynomial / in a; (q^) is invariant under the above involution, it 

is a polynomial in the sinusoidal coordinate ri{x; A): 

X(/») = fix) ^ fix) = f{vix;\)). (2.31) 

The system is shape invariant [^ [2] , 

^(A)^(A)^ = fiAiX + <5)U(A + d)+ £:i(A), (2.32) 

which is a sufficient condition for exact solvability and it provides the explicit formulas for 
the energy eigenvalues and the eigenfunctions, i.e. the generalised Rodrigues formula [2]. 
The forward and backward shift relations are 

TiXjPnix; A) = SniX)Pn-iix; A + (5), S(A)F„_i(a;; \ + d) = P„(x; A), (2.33) 

where the forward and backward shift operators are 

J^(A) = B{0;\)^{x]\)-\l-e^), i3(A) = B{0;X)-\B{x] X)-D{x; X)e-^)^ix; X). (2.34) 

2.3 Symmetry and virtual state vectors 

Let us define the twist operation t of the parameters: 

t(A) = (A4-Ai + 1,A4-A2 + 1,A3,A4), t' = id. (2.35) 



We introduce two functions B'{x) and D'{x) by 

B'{x;X) = B{x;i{X)), D' {x; X) = D {x; i{\)) , (2.36) 



namely, 



[x + d — a + l){x -\- d — b + l){x + c){x + d) 
{2x + d){2x + 1 + d) 



:R 






{(x + a — l)(a; + 6 — l)(x + (i — c)x 
(2a;-l + d)(22; + rf) '" . . 

c^g (1 - ag"-^)(l - 6g"-^)(l - c-^dg")(l - g") ' ^ ^^^^ 

We restrict the parameter range 



: qR 

R 

: gR 



R : d + M <a + b, qR : ah < dq"' , (2.39) 

in which M is a positive integer and later it will be identified with the total number of 
deleted virtual states. It is easy to verify 

B{x; X)D{x + 1; A) = a{XfB'{x; X)D\x + 1; A), (2.40) 

B{x; A) + D{x; A) = a(A) {B\x- A) + D\x] A)) + a (A), (2.41) 

B\x-X)>Q (x = 0,l,...,x^ax + M-l), (2.42) 

D'(x;A)>0 (x = l,2,...,x^ax), I^'(0; A) = D'(x^ax + 1; A) = 0. (2.43) 

Here the constant a(A) is positive and a'{X) is negative: 

/i\ fl ^R « //x\ f — c(a + 6 — d — 1) :R /„,,\ 

< "^^) = I a6rf- V^ :gR' ^ > « ^^^ = { -(1 - c)(l - a^V) : gR " ^2-44) 

The above relations (I2.40p - (l2.43j) imply that we can define a virtual Hamiltonian l-i' by the 
twisted parameters (the A dependence is suppressed for simplicity): 



n'{X)=H{i{X)) = -^/mx) e^ ^U{x) - ^W{x) e-^ ^W{x) + B'{x) + D\x), (2.45) 
and the original Hamiltonian and the virtual Hamiltonian are linearly related 

H(A) = a(A)H(t(A)) + a'(A). (2.46) 



This also means that 'H(t(A)) is positive definite and it has no zero-mode. In other words, 
the two term recurrence relation determining the 'zero- mode' of 'H(t(A)) 



^(t(A))0o(x;A) = O (x = 0,l,...,x„,ax-l), (2.47) 



can be 'solved' from a; = to x = Xmax — 1 to determine 



^o(a;;A) 



Sl^^, (.x = 0,l......».). (2,48) 



But at the end point x = Xmax, the 'zero-mode' equation f l2.47p is not satisfied, because of 
the boundary condition 

5'(x„ax; A) 7^ 0, ^(t(A))0o(a^;A)^O (x = X,nax). (2.49) 

The new Schrodinger equation 

n'iX)Mx; A) = Sl{\)U^; A), (2.50) 

can be almost solved except for the end point x = Xmax by the factorisation ansatz 

4)^{x;X) = (t)oix;\)^^{x;X), 

as in the original (g-)Racah system. By using the explicit form of 4>o{x; A) (12.481) . the new 
Schrodinger equation for x = 0, . . . , Xmax — 1 is rewritten as 

B'{x- A) (ev(x; A) - U^ + 1; A)) + D'{x; A) {Ux; A) - U^ - 1; A)) = £U^)Ux; A). (2.51) 

This is the same form of equation as that for the (g-)Racah polynomials. So its solution for 
X G C is given by the (g-)Racah polynomial (12.261) with the twisted parameters: 

e;(x; A) = A {x; t(A)) , SliX) = S^ (t(A)) . (2.52) 

Among such 'solutions', those with the negative energy and having definite sign 

^v(a;;A)>0 (x = 0, 1, . . . ,Xmax,a;max + 1; v G V), (2.53) 

^,(A) <0 (vG V), (2.54) 



are called the virtual state vectors: {(/)v{x)}, v G V. The index set of the virtual state vectors 



IS 



V = {1,2, 



) ^ max J ) 



}, v^ax = min{[Ai + A2 - A4 - 1]', [|(Ai + A2 - A3 - A4)]}, (2.55) 



where [x] denotes the greatest integer not exceeding x and [x]' denotes the greatest integer 

not equal or exceeding x. We will not use the label state for deletion, see f l3.58p . The 

negative virtual state energy conditions (12.541) is met by v^ax < [^1 + A2 — A4 — 1]'. For the 

positivity of Cv{x; A) fl2.53p . we write down them explicitly: 

■— V, V — a — b + c + d+1, —X, x + d 
d — a + 1, d — b + 1, c 



4-^3 



^v(a;; A) 



4^3 



q ^, a b cdq^'^ , q ^, dq 
a~^dq, b'^dq, c 



q;q 



:R 
: gR 



E(— V, V — a — b + c + d+1, —X, x + d)k 1 „ 
r-. -^—. —. TT :R 



fc=0 



fc=0 



{d — a + l,d — b + l,c)k 
a~^b~^cdq''~^^,q~^,dq^;q)k q^ 



k\ 



(2.56) 



(a ^dq,b ^dq,c]q)i 



(?;?); 



: gR 



Each fc-th term in the sum is non-negative for 2vjnax < Ai + A2 — A3 — A4. 
Here is a summary of the properties of the virtual state vectors: 

(t)Q{x;\) = (j)o(x;i{X)), 0v(a;; A) = (p^{x;i{X)) = if)oix;X)^^{x;X) (v G V), 

^^{x;X) = Pv(x;t(A)), ^v(2;; A) = ^^(^(a;; A); A), 

n{X)4>^{x; A) = £^{X)(f)^{x; A) (x = 0, 1, . . . , x^ax - 1), 

'H(A)0v(a:max; X) ^ £v(A)0v(a;max; A), ^^^(A) = £^{i{X)), 

-{c + v){a + b- d-1 -v) :R 
-(l-cg^)(l -afecf-ig-i-^) : gR ' 

r(l - a)T{x + b)T{d -a + l)T{b - d - x) 

def 00 (2;; X) 



£^{X) = a{X)SUX) + a'{X 
u{x] A) 



(a;;A) 



r(l - a - x)T{b)T{x + d-a + l)T{b - d) 
{a'^q^'-^, b, a'^dq^'^^, bd^^; g)oo 



:R 



(2.57) 
(2.58) 

(2.59) 
(2.60) 

(2.61) 



: gR 



(a ^q,bq'',a ^dq,bd ^q ^■;g)oo 
Note that a'{X) = So{X) < 0. The function z/(x; A) can be analytically continued into a 
meromorphic function of x or g^' through the functional relations: 



, ^ ^, B(x;X) , ., , _, -, D(x;X) , ^ 
v{x + 1; A) = \ L t^{x] A), u{x - 1; A) = ^./ i. v{x; X 



By ^(xmax; A) 



aB'{x]X) 
0, it vanishes for integer x 



(2.62) 



aD'{x; A) 
max + 1 < a; < Xmax + M, u{x; A) = 0, and at 



negative integer points it takes nonzero finite values in general. 
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3 Multi-indexed (g-)Racah Polynomials 

In this section we apply the Crum-Adler method of virtual states deletion to the exactly 
solvable systems whose eigenstates are described by the (g-)Racah polynomials. Since all 
the eigenvalues remain the same, i.e. the process is exactly iso-spectral deformation, the size 
of the Hamiltonian is unchanged. 

Various quantities are neatly expressed in terms of a Casoratian, a discrete counterpart 
of the Wronskian. The Casorati determinant of a set of n functions {fj{x)} is defined by 



W[/i,...,/„](x) = det(/fc(x + j-l)) 



l<j,fc<n 

(for n = 0, we set W[-](x) = 1), which satisfies identities 

n-l 

W[gfi,gf2, ..., gfn]{x) = l[g{x + k)- W[/i, /s, • • • , fn]{x), 

k=0 

W[W[/i,/2,...,/n,^],W[/i,/2,...,/n,/i]](x) 

= W[fi,f2,...,fn]{x + l)W[h,f2,...,fn,g,h]{x) (n>0). 

3.1 Virtual states deletion 

Let us provide the basic formulas starting from one virtual state deletion. For simplicity 
of presentation the parameter (A) dependence of various quantities is suppressed in this 
subsection. 

one virtual state vector deletion 

First we rewrite the original Hamiltonian by introducing potential functions Bdi{x) and 
Ddi{x) determined by one of the virtual state polynomials Cdi{x) (di G V): 

BdAx)'='aB'{x)H^, bdAxV^aD\x)^^. (3.1) 

We have Bd^{x) > (x = 0, 1, . . . , x^a^), i)di(O) = DdA^m^^ + 1) = 0, Dd^{x) > (x = 

1, /, . . . , Xjnax) ^-Ild 

B{x)D{x + 1) = BdAx)DdAx + 1), 
B{x) + D{x) = Bd, (x) + Dd, (x) + Sd, , 

where use is made of (12.511) in the second equation. The original Hamiltonian reads: 

it ^ "^di d,l ' ^d\i 

11 



The virtual state vector (pdii^) is almost annihilated by Ad,, except for the upper end point: 

Ad,4>di{x)=0 (X = 0,l,...,Xmax - 1), ^di<^di(a^max) 7^ 0. 

The proof is straightforward by direct substitution of (13. ip and fl2.48p . 

Next let us define a new Hamiltonian Tid, by changing the order of the two matrices Ali_^ 
and Ad, together with the sets of new eigenvectors (pdini^) and new virtual state vectors 

'Hdi = AdiAj^ + Sdi, T^d, = {T^dxx,y) (3^5 2/ = 0, 1, . . . jXmax)) (3.2) 

4>d,n{.x) = Ad,(t)n{.x) (X = 0,l,...,a;max;'^ = 0,l,...,?2inax), (3.3) 

0div(a;) = Ad,(i)^{x) + (^^.x.^axV'rfiv (X = 0, 1, . . . , Xmax; V G V\{dl}), 



def V^aS'(Xmax) 0o(2;max) c I \C I i lA ^Q/1^ 

V^div = ; =4rfi(3^max)4v(a:max + !)• (3.4) 



y6i(a;max)^di(a;E 



The (/9djv term is necessary for the Casoratian expression for (t)d,w{x) in (13. 5p to hold at 
X = Xmax- It is easy to verify that 0^^,1(2;) is an eigenvector and that (pdiv{x) is a virtual 
state vector 

li-diVdinX,-^ ) ^nVdinX,-^ ) y-^ U, i, . . . , Xjnax; f^ U, i, . . . , JT-maxj i 
'Hd,4>div{x) = S^(f>d,v{x) (X = 0, 1, . . . , Xmax - 1; V G V\{(ii}), 
'T-dir'aiivl^'^inax/ 7^ ^vV^a!ivl^''^niax/- 

For example, 

'Hd,4>din = (Ad^A^^ + ^dj^d,0„ = Ad,{Ali^Adi + £di)(pn 

•^di'l-Vn •^di^nVn ^rf^di^Pn ^n^Pdin- 

The two Hamiltonians "H and Tid, are exactly iso-spectral. If the original system is exactly 
solvable, this new system is also exactly solvable. The orthogonality relation for the new 
eigenvectors is 



'Mm, (Pdim) — 

x=0 



y^ (f}din{x)(l)d,m{x) 
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= {^n-^di){4>n,4>m) = {^n -^dJ-W^nm {n,m = 0, 1, . . . , rimax)- 

This shows clearly that the negative virtual state energy {£^ < 0) is necessary for the posi- 
tivity of the inner products. 

The new eigenvector ^^^^(x) (13. 3p and the virtual state vector (pdivix) (13. 4p are expressed 
neatly in terms of the Casoratian (x = 0, 1, . . . , Xmax) 

0rf,„(x) = -^=^d^^^^w[u,j^Pn]{x), (pdA^) = -^=^d^^^^w[u,mx). 



(3.5) 
We will show that the positivity of the virtual state vector is inherited by the new virtual state 
vector (pdiv{x) (13.41) . The Casoratian W[,^(ij, Cv]{x) has definite sign for x = 0, 1, . . . , Xmax + 1, 
namely all positive or all negative. By using (I2.5ip we have 

aB\x)W[U,^{x) = aD\x)W[U,^ix - 1) + (4 " £.)Uix)Ux). 
By setting x = 0, 1, . . . , Xmax + 1 in turn, we obtain 

±(4 - ^,) > ^ ±w[erf, , e'v] {x)>o (x = 0, 1, . . . , x^a. + 1). 

Note that the set of virtual eigenvalues {4} ^^^ mutually distinct. We will now show that 
the new groundstate eigenvector (t>dio{x) is of definite sign as the original one (poi^) (12. 8p . 
We show that the Casoratian W[^di, H(^) has definite sign for x = 0, 1, . . . , Xmax- By writing 
down the equation 7^0„(x) = £n4>n{x) (x = 0, 1, . . . , Xmax) with T-i = A\ Ad^ +4 ^^^ (I2.14p . 
we have 



aB'{x)u{x + l)P„(x + 1) + aD'{x)iy{x - l)P„(x - 1) = {B{x) + D{x) - £n)iy{x)Pn 



X 



In terms of the functional relations of z/(x) (I2.62p . it is reduced to the original difference 
equation for P„(x) and it is valid any x G C By using this, we can show 

aB'{x)W[U,i^Pn]ix) = aD'(x)W[4, z/P„,](x - 1) + (4 - £„,)4(x)//(x)P„(x). 

By setting n = and x = 0, 1, . . . , Xmax in turn, we obtain 

-^[idi,Jy]{x) >0 (X = 0,1,..., Xmax). 
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Let us rewrite the deformed Hamiltonian T-L^i in the standard form. The potential func- 
tions Bci^{x) and Dd^{x) are introduced: 

Bd,{x) = aB{x+l) — ^, (3.6) 

D,,{x) '^'aD'{x)H^±^'^^^M^^^. (3.7) 

The positivity of B^^ (x) and D^^ (x) is shown above and the boundary conditions B^^ (xmax) = 
and -Drfi(O) = are satisfied. They satisfy the relations 

Ba,{x)Dd,ix + 1) = Bd,{x + l)Da,ix + 1), 
Bd, {x) + Dd, {x) = Bd, {x) + Dd, (x + 1) + £d^ ■ 

The standard form Hamiltonian is obtained: 

nd,=AiAd,, (3.8) 

Ad/=VBdAx)-eyDd,{x), Al = ^Bd,{x)-^Dd,{x)e-^, (3.9) 

in which Ad^ annihilates the groundstate eigenvector 

Ad^<Pdio{x) = (X = 0, l,...,Xmax)- (3.10) 

This one virtual state vector deletion is essentially the same procedure as that developed 
for the exceptional orthogonal polynomials in [23]. See § 13.21 for the explicit expressions. 

multi virtual state vector deletion 

We repeat the above procedure and obtain the modified systems. The number of deleted 
virtual state vectors should be less than or equal |V| and M. 

Let us assume that we have already deleted s virtual state vectors {s > 1), which are 
labeled by {di, . . . , d^} {dj G V : mutually distinct). Namely we have 

'Hdi^As = •^di...ds^di...ds ~^ ^dsi T^di...ds = {'Hdi...dsX,y) (a^,?/ = 0, 1, . . . ,Xmax), (3-11) 



Ad,...d. = \jBd,...dX^) - e'^^Dd^.^dAx), ^^,...,, = \/Bd,...dAx) - \jDd,...dX^) e'^ (3.12) 

6 /Ndef „,, , ,s W[^'d,,...,^'d^_J(x) W[<erf,,...,^"rfJ(x + l) 

Bd^ dAx) = aB [x + s — 1) ; ; ; , 3.13 

bd, d (X) ^=^ aD'(x)M^^^llllM^^±il ^[^\'--^\]^"-^\ (3.14) 

w[e.„...,ed._j(x) w[u,...,a](^) 
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(I)di...dsn{x) = Adi...ds(pdi...d,.in{x) (X = 0,1,. . .,Xmax;"- = 0,1,.. .,nmax), (3.15) 

(i>di...d,A^) = ^di...d>di...d,_iv(a;) +Sa;,x„,^^fdi...dsv, {x = 0, l,...,Xmax;v G V\{di, . . . , c/J), 

^d-,...dsv - (Pdi...dsO[Xmi,x)——~ J TTv^TTTT ? u 7' [O-i-0) 

T-idi...ds(l>di...d,n{x) = Sn(pdi...dsn{x) {x = 0,1, . . . , X^^,^] U = 0, 1, . . . , Tlmax) , (3.17) 

'Hdi...ds4>di...dsvix) = £v(i)dr...dsvix) (X = 0,l,...,Xmax- 1; V G V\{rfi, . . . , 4}) , (3.18) 

(</>di...d.n,0di...d,m) = "^ (pdi...dsnix)<j)di...dsm{x) = JJ(^n-^d,) " -^^nm 

x=0 j=l " 

(n,m = 0,l,...,ninax)- (3.19) 
The eigenvectors and the virtual state vectors have Casoratian expressions (x = 0, 1, . . . , Xmax) 



(-l)%/n;=i«5'(a; + J-l)0o(x)W[e,,,...,?d,,/vF„](x) 
(f>d,...d.n{x) = ^ , , (3.20) 

^JW[U,■■■,L]{x)W[U,■■■,U]{^ + l) 



(-l)%/n;=i«5'(a; + J-l)</'o(2;)W[^,,,...,Cd.,?v](a;) 
0d,...d.v(x) = ^ , . (3.21) 

^Jnu ,■■■, L] (x) w[u , . . . , a] (x + 1) 

The Casoratian in the virtual state vectors W[^d^, . . . ,^ds,^v]{x) has definite sign for x = 
0,1,..., Xmax + 1, and that appearing in the groundstate eigenvector W[^di , ■ ■ ■ ,Cdsi''^] (x) has 
definite sign for x = 0, 1, . . . , Xmax, too. 

The next step begins with rewriting the Hamiltonian 'Hdi...ds by choosing the next vir- 
tual state to be deleted ds+i G V\{di, . . . ,ds}. The potential functions Bd-^^,,,ds+i{x) and 
Ddi...ds+i{x) are defined as in (I3.13j) - (l3.14j) by s -> s + 1. We have Bdi...ds+i{x) > (x = 

0,l,...,Xmax), Ddi...ds+i{0) = i)di...d»+i(a;max + 1) = 0, Dd^...ds+i{x) > (x = 1, 2, . . . , Xmax)- 

These functions satisfy the relations 

Bd^...d,+A^)Dd^...d,+A^ + 1) = Bd,...ds{x + l)Dd^...d.X^ + 1)> 
^di...d,+i(x) + Dd,...ds+Ax) + ^d,+i = Bd,...d.Xx) + Dd,...dAx + ^) + ^ds- 

The Hamiltonian 'Hd^...ds is rewritten as: 

Tidi...ds — •^di...ds+i '^i---"'s+i "T ^ds+l) 



4£.f /r. . (^\_^dlr>, . / ^ ^t 



^di...d,+i = y Bd,...d,^,ix) - e yDd,...d,+,ix), ^d,...d^^, = \^Bd,...d,+Ax) - y Dd,...d,+,ix) e 
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Now let us define a new Hamiltonian l-Ldi...ds+i by changing the orders of Al^_^ ^^ ^ and 
•^di...ds+i together with the eigenvectors (l)di...ds+in{x) and the virtual state vectors (pdL.As+i v{x): 

il-di...ds+i = •^di...ds+1-^di...ds+i ' ^^s+iJ ridi...ds+i — [''-di...ds+ix,y) (,2^) 2/ = U, i, . . . , Xmaxj j 
'Pdi...ds+i n[,-^ ) •^di...ds+i Ydi...ds n[,-^ ) {-^ U,i,..., Xmaxi '^ U, i, . . . , Tlj^axj j 

(pdi...d,+ivix) = Adi...ds+i<l>di...dsA^) + ^x,x,^^^Vdi...d,+iv, 

(X = 0, 1, . . . , Xmax; V G V\{di, ..., ds+l}), 

def , . . W[^d, , . . . , 6,+i] (Xnia^)W[^d, , • • • , 6s , lv] (^max + 1) 
V^di...ds+iv — <Pdi...d,+^0[Xmax) 



W[^rfi , . . . , ^dj (Xmax + l)W[^d, , . . . , ^d^+, , //] (Xmax) 

The orthogonality relation reads 

s+l ^ 

(0di...d,+in,0di...d,+im) = ]^(^n " ^d, ) " "^-'^nm («, "T- = 0,1,.. .,nmax)- 

i=l " 

The functions 0di...ds_j,in(a:;) and (l)di...ds+iv{x) are expressed as Casoratians as in fl3.20p - 
(I3.2ip . The Casoratian W[^di, • • • , ^d^+n Cv]{x) has definite sign 

w[e,,,...,a^j(o)w[e,„...,a,ev](o) 

^ ±W[|rf,, . . .,U+i,iv]{x) > (X = 0, 1, . . . ,Xn,ax + !)• 

Likewise W[^rf^, . . . ,Cds+n H(^) ^^^ ^^^ lowest eigenvector 0di...d^+io(2^) have definite sign 

w[e,„...,a](o) ,Q 

=^ ±W[6,,...,6^+,,Z/](2;) > (X = 0,l,...,Xmax)- 

These establish the s + l case. 

At the end of this subsection we present this deformed Hamiltonian 1-Ldi...ds i^^ the standard 
form, in which the A operator annihilates the groundstate eigenvector: 

'Hdi...ds = '^di...ds-^di...ds, (3.22) 

Ad,...d. = VBd,...dA^) - eWDd,...dA^), A...d. = VBd,...dA^) - ^Dd,...dX^)e-^ (3.23) 

which satisfies 

Ad^...ds4>d^...dsQ{x) = ^ (X = 0,l,...,Xmax)- (3.24) 
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The potential functions Bd^,,,ds{.x) and -Ddi...ds(a^) are: 

n / ^def nV ^ W[g'rf,, . . . , IdJla; + 1) W[U, ■ ■ ■ ,U,'^]{x - I) ,„ .^^ 

-t^di dJa:) = aiJ (a;) ; ; ; . [6.2b) 

The positivity of Bd^,,,d,ix) and Dd^...ds{.x) is shown above and the boundary conditions 
-Bdi...ds(a^max) = and Dd^...dsiS^) = are satisfied. They satisfy the relations 

5di...d.(a;)I^di...d,(a; + l) = Bd^...ds{x + l)Dd^...ds{x + 1), 
Bd^...ds{x) + Dd^...dX^) = Bdi...dA^) + Dd^...dA^ + 1) +^d,- 

It should be stressed that the above results after s-deletions are independent of the orders 
of deletions {(pdi...dsn{x) and (t)di...dsv{x) may change sign). 

3.2 Explicit forms of mult i- indexed (g'-)Racah polynomials 

Here we present the main results of the paper. The multi- indexed (q'-)Racah polynomials 
are obtained by applying the method of virtual states deletion to the (g-)Racah system. 
The parameter A = (Ai,A2,...) dependence is now shown explicitly. The eigenvectors of 
the models in § 5 of [2\ are described by orthogonal polynomials in the sinusoidal coordinate 
ri{x; A). The auxiliary function ip{x; A) is defined by 

.^,rKx+2|A)_rKx^^ (3.27) 

r/(l;A) 



and it satisfies (with d defined in (I2.15p - (l2.16p ) 

ip{x; A) 



(^(1;A). (3.28) 



ip{x - 1; X + 26) 
All the models in § 5 of [2] have shape invariance [44j- The following relations are very useful: 



B(0;A) M^-A + S) g(0;A) 0o(x; A + (5) 

^^'-^^^ = ^lBi^) 0o(x;A) ' ^^^^^^ = iwTT^)M- + i;><y ^'-''^ 

B{x + l;X) '^ <^(x;A) ' ^(a;; A) '^ ^{x\X) ' ^ ' ' 

We delete M virtual state vectors labeled by 

V = {di, d2, ■ ■ ■ , rfjv/} {dj G V : mutually distinct), (3.31) 
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and denote Hdi...dM, 4'di...dMn, Adi...dM, etc. by Kv, 4>vn, Av, etc. 

Let us denote the eigenvector (f)vn{x; A) in f l3.20p after M deletions (s = M) by 0|^^(x; A). 
We define two polynomials Hx)(x; A) and -Pp,„(x; A), to be called the denominator polynomial 
and the multi-indexed orthogonal polynomial, respectively, from the Casoratians as follows: 



W[^di, . . . , CdM]i^i ^) = Cv{X)Vm{x\ A)S©(x; A), 

W[^d,, . . . , idM,^Pn]{x] A) = Cv,nW'^M+l{x] X)Pv,n{x] \)p{x] A + M~5), 



7. def 



(5= (0,0,1,1), t(A) + /?5 = t(A + /?(5) (V/?G 



(3.32) 
(3.33) 
(3.34) 



The constants C-p(A) and Cx),„(A) are specified later. The auxiliary function lpm{.x;X) is 
defined by |12|: 

rj{x + A; — 1; A) — rj{x + j — 1; A) 



lPm{.X]\) = Yl 



v{k - j; A) 



l<j<k<M 

and i^oix; A) = ipi{x; A) = 1. The eigenvector (I3.20p is rewritten as 



(3.35) 



^r„(^; A) 



^\M \M(M-l) ^T^,nW 



X 



bo{x- A + MS) 



At 



l[a{X)B'{0;X+U-l)~S) 

Pv,n\X', A). 



(3.36) 



'Ev{x;X)Ev{x + 1; A) 
The multi- indexed orthogonal polynomial Px>,n{x] A) (13.331) has an expression 



Pv,n{x] A) = Cr.,„(A) ^(pM+i{x] A) ^ 

(di(2;i) ■■■ {dMi.xi) 

idAx2) ■■■ ldM{x2) 



X 



ri{xi)Pn{xi) 

r2{x2)Pn{x2) 



(3.37) 



■^^1(2^*/+!) ■■■ idM{xM+l) rM+l{xM+l)Pn{xM+l) 

where Xj =' a: + j — 1 and rj{x) = rj{x] A, M) (1 < j < M + 1) are given by 

{x + a,x + h)j_i{x + d- a + i,x + d-h + j)m+i-j 



def 



rj{x + 3-l;X,M) 



def J {d — a + l,d — b + 1)m 



:R 



{abd-^q-^y-^q^^''{a-^dq, b-^dq; q) 



M 



: qR 



(3.38) 
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M 
def ■ 



One can show that E-d (I3.32p and Pv,n (I3.37P are indeed polynomials in rj: 

Ev{x; A) = Ej,{r]{x; A + (M - 1)^); A), P^,„(x; A) = Pv,n{v{x; A + M^); A), (3.39) 
and their degrees are generically £ and i + n, respectively (see (I3.59P ). Here i is 

M 

J2dj-lM{M-l). (3.40) 

i=i 

The involution properties (I2.3ip of these polynomials are the consequence of those of the 
basic polynomials -P„(x) and Cdj{x). We adopt the standard normalisation for Ex> and Pv,n'- 
S-p(0; A) = 1, -pD,n(0; A) = 1, which determine the constants ^©(A) and Cx),n(A), 

r rv^42.f 1 TT gd,(A)-grf,(A) 

y?Af+i(0; A) ^^-^ a[X)B'[j - 1; A) 

The use of dual polynomials Qx{£n) = Pn{v{^)) [2] is essential for the derivation of these 
results. The three term recurrence relations of {Qx{^)} are specified by B{x) and D{x). The 
denominator polynomial Hx)(a;; A) is positive for x = 0, 1, . . . jXmax + 1- The lowest degree 
multi-indexed orthogonal polynomial Pv,o{x; A) is related to Sx)(x; A) by the parameter shift 
X^X + d: 

Pv,o{^-X) = Ej,{x;X + S). (3.43) 

The potential functions i?p and Dx) (I3.25|) - (I3.26|) after M-deletion [s = M) can be expressed 
neatly in terms of the denominator polynomial: 

g.(x; A) ^ B(x; A + MS) ^ Y^' f ' -"'^ ^- ^+/' . (3.44) 

D.(.; A) ^ D(.; A + Mg) "^'^ \^'> 't 7 \ ^ t'^ ■ (3.45) 

These formulas look similar to those in the exceptional polynomials [23]. The groundstate 
eigenvector (/)x>o is expressed by (l)o{x) (12. 8p and Hx)(2;; A): 



y© 



o(a;; A) 



\y=o^»(?/ + l) V ^p(2;;A)^i,(a; + l;A) 
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= iIjv{x; X)Pv,o{x; A) oc (;/)|,™(x; A), (3.46) 

Mx;^) = \/^A^) : ^°^"'^ + '^^^ , V'.(0;A) = 1. (3.47) 



'Sx,(x; A) Hx)(x + 1; A) 
We arrive at the normalised eigenvector (j)vn{,x; A) with the orthogonahty relation, 

(t>vn{x\ A) = tpvix] \)Pv,n{x\ A) (X 0|,^^(a;; A), 0d„(O; A) = 1, (3.48) 

V ^^ , ' ^. Pv,n{x;X)Pv,m{x-\) = J" (n, m = 0, 1, . . . , n^ax) • (3.49) 

^ ^d(1;A) (i„(A)^(i75,„(A)^ 

It is worthwhile to emphasise that the above orthogonality relation is a rational equation of 
A or g"*^, and it is valid for any value of A (except for the zeros of denominators) but the 
weight function may not be positive definite. 

The shape invariance of the original system is inherited by the deformed systems. The 
matrix Adi...ds+i{^) intertwines the two Hamiltonians T-Ldi...dsW and T-Ldi...ds+iWy 



Adi...ds+iW -Adi^As+iW - 'Hdi...dsW -£ds+iW, 

Adi...d,+iWAdj_...ds+iW = ^di...ds+iW - ^d,+i(A), 

and it has the inverse. By the same argument given in §4 of [19j, the shape invariance 
of 'H(A) is inherited by 'Hd^(A), 'Hdid2W, " " " • Therefore the Hamiltonian ^^^^(A) is shape 
invariant: 

AvWAviXy = ^^Avi^ + 5)Up(A + S)+ £i(A). (3.50) 

As a consequence of the shape invariance and the normalisation, the actions of Axi{X) and 
AT){\y on the eigenvectors (/)T:>n{x; A) are 

£ (X) 

AvW4>Vn{x; A) = — ^ " (pVn-l{x; X + S) (X = 0, l,...,Xmax- 1), (3.51) 

^-^^^(O; A) 
AT){Xy(pVn-iix; X + d) = a/ 5^,(0; A) (pvnix; A) (x = 0, 1, . . . , x^ax)- (3.52) 

The forward and backward shift operators are defined by 



JS(A) = VBv{0; A) ilJv (x; X + 6)-' o ^^(A) o ^^ (x; A) 



Lp{x; X + MS)'^T>{x + I; X' 



■ (Evix + l;X + 6)- Ev{x; X + <5)e^) , (3.53) 



def 1 



Bv{X) "=^* , ^ ^^ ^v (x; X)-' o Av{Xy o ^^ (x; X + 6) 
y/Bx){0; A) 
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= ^-. (3.54) 

X (b{x; X + MS)^T)ix; X) - D{x- X + M~5)'^v{x + 1; A)e-^V(x; A + M~5), 

and their actions on Px),n(a^; A) are 

J'v{X)Pv,n{x-, A) = £n{X)Pv,n-i{x-, A + (5), Bv{X) Pv ,n-i{x-, X + 5) = Pv,n{x-, A). (3.55) 

As in the original (g-)Racah theory (I2.33p . these formulas are useful for the explicit calcula- 
tion of the multi-indexed polynomials. The similarity transformed Hamiltonian is 

nv{X) = ^T,{x- A)-i o UviX) o ^^{x- A) = Bv{X)Tt,{X) 

Hx)(x;A) /^Sx)(x + 1; A + (5) ^q 



= B{x; X + MS) , "^ ' ^ .' ' ' - e' 

Ev{x + l;X)\ Sx)(x;A + <5) 

+ C(x;A + Ma)i#^±ii^(i#i=^iA±i)-e-4 (3,56) 

^T>{x;X) \ r.T>{x;X + d) ) 

and the multi-indexed orthogonal polynomials Pv,n{x] A) are its eigenpolynomials: 

•H©(A)A?,„(x; A) = £n{X)PvA^; A). (3.57) 

Other intertwining relations are 

n'^Ad,...d,+,iX + S)Ad,...dsW = ^di...d,+i(A)^di...d,+i(A), 
K"Mdi...d,+i(A)^di...d,(A)"f = ^di...d,+i(A)Ud,...d,+,(A + (5), 

with the potential functions given in (I3.13p - (l3.14p (with s — )■ s + 1). 

Including the level deletion corresponds to M — 1 virtual states deletion: 

A?,n(a:;A) =Pv',n{^;X+~d), V = {d^- 1,. . . ,dM-i- !}■ (3.58) 

dM=0 

This formula is similar to those in the multi-indexed Jacobi theory, eqs.(48)-(49) in [T]. The 
denominator polynomial S© behaves similarly. This is why we have restricted dj > 1. 
The coefficients of the highest degree terms of the polynomials Sx) and Pv,n, 

Sx)(2/; A) = c^{X)y^ + (lower order terms), 
Pv,n{y] A) = Cx>^ni^)y ^" + (lower order terms). 
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are 



jM 



c|(A) 



n =i(-« -h + c + d + dj + l). 



M 



Hi 



<j<k<M 



'—a — b + c + d + dj + dk + 



jtUcT^-^cdq'^^+^q),^ 



^ i=l 



(c, (i — a + 1, (i — 6 + l)j-i 
{c,d — a + l,d — b + l)d 



M 



-^U 



{c,a ^dq,b ^dq]q)j_i 



ni<,<fe<Af (1 - a-^b-^cdq<^^+<^^+^) AA (c, a-^dg, ft-irfg; g), 
{a + b + c- d + n-l)n{c)M 



iM 



C©,n(A) = C|(A) X < 



(a, &,c)„nj=i(c + "' + c?j; 
{abcd~^q'^~^;q)n{c;q)M 



R 
gR 



:R 



: qR 



(3.59) 



The exceptional X^ (g-)Racah orthogonal polynomials presented in ^23j correspond to 
the simplest case M = 1, V = {£}, i > I: 



^e{x; A) = H{£}(x; A + M - 6), Pe,n{x; A) = P{£}^n{x; X + id - 6). 



(3.60) 



4 Summary and Comments 

Following the examples of multi-indexed Laguerre and Jacobi polynomials [1], multi-indexed 
(g-)Racah polynomials, the discrete quantum mechanics counterparts, are constructed. These 
new polynomials could be considered as a further generalisation of Bannai-Ito polynomials 
[15] . The next stage will be the construction of multi-indexed Askey- Wilson and Wilson 
polynomials. The basic logic is the same for the ordinary quantum mechanics as well as for 
the discrete quantum mechanics with real [2] or pure imaginary shifts [28]. Starting from 
the factorised Hamiltonians of exactly solvable quantum mechanical systems, a series of new 
'deformed' exactly solvable quantum systems are generated by applying Crum-Krein-Adler 
formulas [36| [37] or multiple Darboux transformations [39] through deletion of various vir- 
tual states instead of eigenstates. The virtual state vectors are polynomial 'solutions' of 
a virtual Hamiltonian which is obtained by twisting the discrete symmetry of the original 
Hamiltonian. They fail to satisfy the Schrodinger equation of the virtual Hamiltonian at one 
of the boundaries, at x = Xmax- When there is only one extra index V = {£} {i > 1), the 
multi-indexed (g-)Racah polynomials reduce to the exceptional polynomials [23l [26]. Like 
the exceptional polynomials, the multi-indexed (g-)Racah polynomials do not satisfy the 
three term recurrence relations. On the other hand, their dual polynomials satisfy the three 
term recurrence relations because of the tri-diagonal form of the Hamiltonian. As for the 
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parameter ranges in which (I2.42p -( l2.43p . f l2.53p -( l2.54p are satisfied, we have taken conser- 
vative ones, f l2.18p -f l2.19p . f l2.39p . fl2.55p . It is quite possible that the valid parameter ranges 
could be enlarged. The difference equations for the multi-indexed (g-)Racah polynomials, 
f l3.55p . fl3.57p are purely algebraic and they hold for any parameter ranges. 

As in the ordinary Sturm-Liouville case (the oscillation theorem) the multi-indexed or- 
thogonal polynomial P-o^niy] A) has n zeros in the orthogonahty range, < y < ?7(xmax; A + 
M6). This is a general property of the eigenvectors of a Jacobi matrix of the form (12. ip . 
fl2:45|) etc. See [16]. 

A few words on the mirror reflection with respect to the mid point of the x-grid, x — )■ 
Xmax — X = N — X. The original (g-)Racah polynomial under the reflection is described by 
the same polynomial with different parameters: 

Pn{N - x; A) = A X P„ (x; (Ai, Ai + A3 - A4, Ai + A2 - A4, 2Ai - A4)) , 

{{a + h - d,a + c- d)n _ „ 

/ _i An (a6rf"\ac(i~i;g)„ 
(" '*' (t.c;,)„ -"^ 

This corresponds to the mirror reflection formula of the Jacobi polynomials: 

Pt'^\-V) = {-irPi^'''\ri), r/(x) = cos2x, n{l - x) = -ri{x) . 

The type 1 and 11 virtual state wavefunctions for the Jacobi case have a twisted boundary 
condition at x = | and x = 0, respectively [1] and their polynomial parts are related by 
this mirror reflection [H]. The virtual state vectors, which fail to satisfy the equation at 
X = ^xnaxi correspond to the type 1. By the mirror reflection, one can consider the type 11 
virtual state vectors, which fail to satisfy the equation at x = 0, instead of x = Xmax- By 
using these type 11 virtual state vectors, the mirror reflexed version of the multi-indexed 
(g-)Racah polynomials can be constructed. They are related to the multi-indexed (g-)Racah 
polynomials by 

F™(x; A) = A X Fi,,„ [N - x; (Ai, Ai + A3 - A4, Ai + A^ - A4, 2Ai - A4)) , 

A-^ = Pv,n{N; (Ai, Ai + A3 - A4, Ai + A2 - A4,2Ai - A4)). (4.2) 

The normalisation conditions are F„(0; A) = Pv,niO; A) = P™"°''(0; A) = 1. Contrary to the 
Jacobi case [1], the type 1 and 11 virtual state vectors cannot be used together to generate 
new multi-indexed (g-)Racah polynomials. 
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Various orthogonal polynomials are obtained from the (g-)Racah polynomials in certain 
limits. Similarly, from the multi-indexed (g-)Racah polynomials presented in the previous 
section, we can obtain the multi-indexed version of various orthogonal polynomials, such 
as the (g-)Hahn, dual (g-)Hahn, alternative g-Hahn (see § V.C.I of [2jj), etc. The infinite 
dimensional cases, the little g-Jacobi, (g-)Meixner, etc will be reported in a separate publi- 
cation. In certain limiting processes, the mirror reflection does not commute with the limit, 
as is well known in the Jacobi — )■ Laguerre limits [HI [1] . The mirror refiexed multi- indexed 
poljTiomials are supposed to play certain roles in such limits. 

Let us emphasise that the discrete symmetries of the original (g-)Racah systems and their 
twisting, which are essential for the construction of virtual Hamiltonians and virtual state 
vectors, are easily recognised in the present parametrisation fl2.20p - fl2.2ip [2], but rather 
unclear in the original parametrisation [H HJ [5]. This is a good reason to promote the 
(g-)Racah systems in our parametrisation. 
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